Non-adiabatic pumping of single electrons affected by magnetic fields 
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Abstract 

On 

Non-adiabatic pumping of discrete charges, reaUzed by a dynamical quantum dot in an AlGaAs/GaAs heterostructure, is studied 
under influence of a perpendicular magnetic field. Application of an oscillating voltage in the GHz-range to one of two top gates, 
crossing a narrow wire and confining a quantum dot, leads to quantized pumped current plateaus in the gate characteristics. The 
Qh regime of pumping one single electron is traced back to the diverse tunneling processes into and out-of the dot. Extending the 
theory to multiple electrons allows to investigate conveniently the pumping characteristics in an applied magnetic field. In this way. 



in 



a qualitatively different behavior between pumping even or odd numbers of electrons is extracted. 
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, 1. Introduction 

Well defined pumping of discrete charges using an excitation 
frequency / has attracted much interest in an eff'ort to realize 
a quantum standard for the ampere dO. Amongst others, in- 
tegrated single-electron circuits (3] as well as quantum infor- 
mation processing |4 | have been proposed. Many approaches 
to obtain such a quantized current have been tested, following 
the pioneering work two decades ago E [6l [71 . In all attempts 
the acid test for a well defined current is the current plateau - 
as a function of some controlling voltage applied at a gate - at 
current values of / = nef. n is the integer number of electrons 
with elementary charge e pumped at frequency /. 

Combining the need for a high current value and high accu- 
racy remains a big challenge. Therefore, devices suitable for 
parallelization and hence multiplication of the current output 
have been put forward in two recent approaches. Beside the 
turnstile device introduced by Pekola et al. |8 1, implemented in 
metallic systems with fixed tunnel barriers, another approach 
uses tunable tunnel barriers in a semiconductor nanowire 191 , 
operated by a single voltage parameter fT0l[TT|[T2l[T3l. In re- 
cent works on AlGaAs based devices of the latter type it was 
shown that a perpendicular applied magnetic field significantly 
modifies the region where quantization appears as well as the 
slope of the quanzited plateaus 114^,15]. The origin of this eff'ect 
is still not understood, but could potentially be used to improve 
the accuracy be several orders of magnitude fTT] . 

Here we off'er a first step towards quantification of this de- 
pendence by applying a recently proposed model for dynam- 
ical quantum dot (QD) initialization |16|. This model links 
the plateaus shape to decay rate ratios between diff'erent charg- 
ing states of the QD and therefore connects it to microscopic 
properties. To this end we will first analyze a relevant set of rate 
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equations for the complete parameter range and then find a re- 
gion in which the above model is valid. In the second part, 
the model is used to extract an even-odd characteristic of the 
plateaus as a perpendicular magnetic field is applied. 

2. Experimental Details 

A picture of one of the two devices used for our measure- 
ments is shown in Fig. [ij^)- Three metallic top gates, 100 
nm in width, are crossing a wire etched in a ^-type AlGaAs 
heterostructure. The distance between the gate centers is 250 
nm and the bottom gate is not used and grounded. Two de- 
vices (labeled 1 and 2) diff'er in channel width, wi = 800 nm, 
W2 = 900 nm, and carrier density Ug^i = 2.1 x 10^^ cm"^ and 

2 = 2.8 X 10^^ cm"^. All measurements were performed in 
a ^He/^He dilution cryostat at temperatures below 50 mK. For 
quantized charge pumping the devices have to be operated in 
the following way: Applying sufficiently large negative volt- 
ages to gate L and R forms a QD between them. An additional 
sinusodial modulation of frequency / is then coupled to gate L. 
If the amplitude is high enough, the energy of the lowest quasi- 
bound state 6o in the QD is pushed below the Fermi energy of 
the leads £'f in the first half of the cycle and can be loaded with 
electrons from source. During the second half of the cycle the 
left barrier is raised fast enough to prevent the electrons from 
tunneling back to source. At the same time eq is raised above 
£^F, so the electrons can be unloaded to drain. This results in a 
current created without an external bias voltage as sketched in 
FigfTJb). 



3. Quantized pumping regime 

To realize transport of a single electron during one oscillation 
cycle, it is necessary to keep the tunnel couplings of the barriers 
in a regime, where exactly one electron is loaded into the QD 
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Figure 1: (Color online) SEM picture of the device shown in (a). Gate voltages 
are indicated, showing gate L colored in red as being modulated. (S) and (D) 
mark soure and drain. The hatched regions are depleted of the 2D electron gas, 
defining a wire. A quasi-bound state is formed between gates L and R, indicated 
by the white ellipse. The third gate is not used and grounded, (b) Schematic 
of the potential along the channel during loading (left) and unloading (right) of 
the quasibound state ^(x). 



from source and subsequently unloaded to drain. First we will 
define a quantization region were this condition is fulfilled 1 17 1. 

The voltages creating the barriers L and R are represented 
by dimensionless parameters Vl and Vr, respectively. Vl(0 = 
+ y^"^^ cos(6i;0 is composed of a DC off'set voltage and a 
radio frequency signal with amplitude y^"^^. It can be assumed 
that the tunnel couplings Y\Vl, Vr) and T'Cy^, Vr) of the left 
and right barrier respectively depend exponentially on the clos- 
ing potential barrier height 1 18 1. In Ref. 1 10| it is shown that the 
resulting time dependence, not only of log but also of log 
can be well approximated by a cosine function, taking into ac- 
count that tunneling through one barrier is additionally aff'ected 
by the other barrier. Hence we assume in dimensionless units: 



Y\yL,VR) = t^V[Vdt)-yVRl 

nyL,y/?) = exp[y;,-ry^(0], 



(1) 

(2) 



where < y < 1 is the cross talk ratio. The value of 6o is 
determined by both voltages in a symmetric way: 



6o(0 = -(ydt) + Vr). 



(3) 



Without loss of generality we can set / = \ m cd = Inf. The 
relative position of the 6o with respect to the Fermi energy is 
fixed by the condition 



60 + ^0 = 0. 



(4) 



Here is a dimensionless parameter which represents the 
Fermi energy in the leads. Thus the shape of the quantization 
steps is completely fixed by the parameters y, and y^"^^. Re- 
garding the temperature T the following assumptions are made: 

i) T is much smaller than the energy modulation amplitude. 
This allows us to replace the full Fermi function by a step 
function. 

ii) T is much larger than the tunneling rates and the modula- 
tion frequency. This allows us to use a simple rate equa- 
tion. IlIOJ 



The level crossing condition ^ is satisfied in the time moments 
t = to -\-m and t = -to -\- m, where m is an integer and 



1 ^0 - - y^^ 
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to = — arccos ^,^p 



(5) 



The range of (V^^, Vr) at which to is real (for fixed V^^) de- 
termines the region of parameters where the current generation 
is possible within the model. During one period from t = -to 
till t = I - to the energy level is below the Fermi surface for 
-^0 < ^ < +^0 and above it for -h^o < ^ < 1 - ^o- Corresponding 
rate equations can be formulated as: 



= {Y^^r)[l-P(t)l -to<t<^to, 



dP 

tunneling into the QD is possible, out is prohibited. 
dP 

— = -(r^ + r'^) [P(0],for + ^o<^< 1-% 

tunneling into the QD is prohibited, out is possible. 



(6) 



(7) 



Imposing a periodic boundary condition in time makes the 
periodic stationary solution P{t) unique. The instantaneous cur- 
rent flowing into the QD from the lead a = /, r is calculated as 



r^(0 dP{t) 



P(i) -h r^(o dt 



(8) 



Corresponding to this model current measurements in lead a 
would ideally give 



= [ F 
Jo 



(0 dt 



(9) 



We will now discuss the parameter range of Vr and V^^ 
where efficient transportation of electrons is possible, shown 
in Fig. |2ja) ifTTll . As sketched by the four pictograms in the 
corners, the heights of the left and right barrier decrease with 
increasing the corresponding applied voltage to the top gates. 
The blue line L indicates F^ = 1 i.e. the transition between the 
left barrier being transparent (above, F^ » 1) or opaque (below 
F^ <^ 1) for charge transfer within one period 1// . In a sim- 
ilar manner, the green line indicates F^ = 1 i.e. the transition 
of the right barrier from transparent (below, F^ » 1) to opaque 
(above, F^ ^ 1). The tilt of the lines relative to the axis is de- 
termined by the cross talk factor y ^ 0. While the lines L and 
R are related to the tunnel couplings the red line D is related 
to the position of 6o relative to the Fermi energy, see Eq. |4] At 
line D 6o crosses the Fermi energy (eo = -^o)- Above this line 
the level is below the Fermi energy and thus the QD is loaded if 
tunnel coupling to at least one lead is sufficient large (F^ » 1). 
In contrast, below line D it becomes energetically favourable to 
unload the QD if it is occupied and sufficiently coupled. 

Note that until now, we have considered only the DC part of 
the applied voltage. In Fig. |2] the rf-modulation of the left bar- 
rier is taken into account by a parallel vertical shift of the three 
lines in both directions. The magnitude of the shift is V^^, as 
indicated by an arrow on the right. The two vertically shifted 
D-lines are marked D+ and D_ . Let us first consider the loading 
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Figure 2: (Color online) (a) The region of Vr and V^^ where efficient trans- 
portation of electrons from source to drain is possible is shaded in gray. The 
derivation and explanation are given in the text, (b) Typical characteristic 
of pumped current of device 1 measured as a function of V^^ and Vr at 
/ = 500 MHz, P = -25 dBm and P = -23 dBm. While the lower and up- 
per line shift with increasing power, the left and right line stay fixed. 

of the QD from the source. Due to the additional AC compo- 
nent loading the QD is possible even for V^^ and Vr voltages 
down to D_ during one part of the cycle. For quantized current 
pumping the loading must occur from the source i.e. via the left 
barrier. Therefore the crossing point of the D_ line and the L_ 
line marks one extremal point (d) of the pumping region. Only 
above D_ and above L_ loading from the source is possible. 
The upper boundary of the pumping parameter region can be 
found when considering unloading the QD to the drain. There- 
fore the line D+ must be considered. Only below D+ unloading 
of the QD is energetically allowed. Additionally for unloading 
to the right drain the right barrier has to be transparent. Thus 
the other extremal point of the pumping regime is marked by 
(b) i.e. the crossing point of D+ and R+. Again, only below D+ 
and R+ unloading to the drain is possible. Connecting the ex- 
tremal points and the lines L_, R+ by vertical lines defines the 
AC pumping region (abed). 

Qualitatively, this region is bordered as follows: At step B the 
QD has been completely loaded exclusively from source during 
-^0 < ^ < ^0- The current outside the pumping region drops 
because the QD unloads to source again instead to drain during 
^0 < ^ < 1 - ^. Outside step R+ the current drops because of 
unsufficient unloading to drain during < ^ < 1 - ^o- Outside 
step F loading during -^o < ^ < takes place predominantly 
from the drain instead from source. Since the electron is also 
emitted to drain during < ^ < 1 - the average dc current 
drops to zero. Outside step L_ the QD is unsufficiently loaded 
through the source barrier, leading again to a current drop. 

In Fig. [2jb) measurements of device 1 are plotted for dif- 
ferent rf-amplitudes. The size of the pumping region (marked 
by dotted lines) increases along the V^^-axis with increasing 
rf-amplitude, but the position of the left and right line remain 



fixed. This eff'ect agrees with the experimental results, hence 
all four lines B, F, L_ and R+ can be identified. The different 
slope at the L_ step might result from neglecting excited states. 
Lowering the right barrier and consequently eq may open ad- 
ditional transport channels through the left barrier into excited 
states. In contrast to the experimental data where several quan- 
tized plateaus appear, only one electron is considered in this 
model. To describe the multi-electron case, for simplicity only 
the transition region across line B will be considered next. 

The transition regions of L_ and R+ can in principle be arbi- 
trarily well separated from those of B and F by choosing a suf- 
ficiently large modulation amplitude. In particular for the mea- 
surements considered below, a range of V^^ has been identified, 
where the transition across line B is not significantly aff'ected by 
that of L, R and F. Therefore, the relevant processes for devi- 
ation from the quantized current value are tunneling events of 
the previously loaded electron back to source. Whenever the 
electron remains on the QD it will be emitted to drain during 
to < t < 1 - so we only have to focus on -^o < ^ < ^o- At 
this time F'' can be assumed. Hence, we can simplify our 
problem to a QD, which decays to source only and becomes in- 
creasingly isolated. This picture can be generalized to describe 
more than one electron, such that the final charge state deter- 
mines the average number of electrons being transfered. 

Such a dynamical QD has already been theoretically investi- 
gated in Ref . |[T6ll : A decay cascade model is proposed in which 
the confining potential is characterized by a fixed sequence of 
decay rate ratios r^/r|^_^ = e^", one rate per each QD occu- 
pation number n. The rates F^ are assumed to decrease uni- 
formly and exponentially in time t and control gate voltage, cor- 
responding to Vr above. This leads to an analytic solution for 
the occupation-number probability distribution which is param- 
eterized by a set of ^„'s. An explicit fitting formula for the aver- 
age number of captured electrons has been derived for ^„ » 1 : 

(n> = |;.--"'"^^' where At = j]6, (10) 

Thus every plateau is described in length and slope by a signif- 
icant S, whereas Sn relates to the (n - l)th plateau. Note that 
apart from the dimension a is identical to the cross talk ratio y. 

4. Experimental results 

Changes in the current plateaus caused by a perpendicular 
magnetic field have already been reported (141 [IS- From the 
investigations above we are able to visualize the changes sepa- 
rately for every plateau as a function of magnetic field. 

In Fig. [3] measurement results are shown. The voltage V^^ 
has been set to the middle of the plateau so that the vertical 
transition regions (corresponding to lines R+ and L_) have no 
siginificant effect on the current, i.e. insufficient filling or emp- 
tying of the QD is avoided. Due to the different carrier densities 
of the two devices, the required power for quantized pumping 
varies. For device 1 the power was set to P = -16 dBm and 
for device 2 to P = -11.5 dBm. The graph shows Sn obtained 
by fitting the normalized current to Eq. [TO] as a function of the 



(^) Device 1 , p>Mr^'^ 


(b) Device 2 

/ 


rr::;;; ^^^^^ 



°0,0 0,5 1,0 1,5 2,0 2,5 0,0 0,5 1,0 1,5 2,0 2,5 



BinT BinT 

Figure 3: (Color online) 6n as function of magnetic field strength of two dif- 
ferent devices. Device 1 (a) was operated at a pumping frequency / = 50 MHz 
and a power of P = -16 dBm. Device 2 (b) was operated at a / = 100 MHz 
and P = -1 1.5 dBm, (b) is shifted for clarity. 

magnetic field. Despite the diff'erence in the wafer character- 
istics, both devices show a reproducible even-odd behavior of 
the slopes. For even n the corresponding ^„'s grow significantly 
faster with magnetic field than for odd n. 

Device 1 shows for the first plateau an increase from d2{B = 
0) = 8.6 to S2(B = 2.75T) = 13.2. This means an increase 
of the tunnel coupling ratio F^/r^ ^ by a factor of circa 100. 
Examining the third plateau the increase of F^/F|^_^ is not as 
pronounced, but still raised by a factor of 15, due to the increase 
from 54(B = 0) = 3.4 to 64(8 = 2.75T) = 8.6. The second 
and fourth plateau show only a marginal increase. The tunnel 
coupling ratio is only raised by a factor of about 2.3 and 3, 
respectively. Device 2 reveals a similar result. While F^/F|^_^ 
increases with magnetic field by a factor of 200 for the first and 
65 for the third plateau, for the second and fourth plateau there 
is no significant increase visible. According to Ref. [1I6II the 
increase of the first plateau of both devices corresponds to an 
enhancement of the accuracy by two orders of magnitude. 

There are several factors contributing to the observed behav- 
ior. Energetic shifts should aff'ect the decay rate ratios. Within 
a single-particle mean-field approach, an even-odd asymmetry 
in the field dependence of decay rate ratios may indicate se- 
quential filling of orbitals. Ellipticity of the confining potential 
as well as deviations from constant-interaction approximation 
may complicate the B -field dependence of orbital energies |[T9ll . 
Also the magnetic field dependent barrier shape has to be kept 
in mind, due to its influence on screening. As a next step a 
microscopic theory of the tunnel couplings has to be found in 
order to obtain an exact information about the origin of this 
even-odd dependent evolution. This insight should be valuable 
to find strategies for systematic accuracy improvements for this 
pumping scheme. 
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